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Abstract. 

We investigate, within Floquet theory, topological phases in the out-of-equilibrium system 
that consists of fermions in a circularly shaken honeycomb optical lattice. We concentrate on 
the intermediate regime, in which the shaking frequency is of the same order of magnitude as 
the band width, such that adjacent Floquet bands start to overlap, creating a hierarchy of band 
inversions. It is shown that two-phonon resonances provide a topological phase that can be 
described within the Bernevig-Hughes-Zhang model of HgTe quantum wells. This allows for 
an understanding of out-of-equilibrium topological phases in terms of simple band inversions, 
similar to equilibrium systems. 
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The recent realisation of topoplogical states in condensed-matter systems (Konig et al. 2007, 
Hasan & Kane 2010, Qi & Zhang 2010) has sparked a flurry of activities in the field of 
cold atoms, aiming at reproducing, engineering, and manipulating these fascinating quantum 
states in traps (Roncaglia et al. 2011, Corman et al. 2014, Burrello et al. 2015) and in optical 
lattices (Hemmerich & Morais Smith 2007, Aidelsburger et al. 2013, Miyake et al. 2013, Jotzu 
et al. 2014, Aidelsburger et al. 2015, Chomaz et al. 2015). 

Topological fermionic systems exhibit protected metallic states at the boundary, while the 
bulk of the material remains insulating. In cold atoms, quantised transverse density currents 
play the role of the quantised charge currents, since we deal with neutral atoms in selected 
hyperfine states, instead of spin-full electrons. Similarly, the pair of atomic hyperfine states 
yields a spin-1/2 structure that allows for the analogue of quantised spin currents obtained 
in condensed-matter systems. The conductivity at the boundary is given by a topological 
invariant, which is quantised and stable against perturbations. 

There are a variety of topological states known by now, and they are well described by 
the ten-fold classification (Ryu et al. 2010), which defines a topological invariant according 
to the symmetries and dimensionality of the system. Although well established, the ten-fold 
classification neglects several effects. First, it does not take into account interactions, which 
leads to the fractional quantum Hall effect (Tsui et al. 1982, Laughlin 1983), to the quantum 
anomalous Hall effect (Nandkishore & Levitov 2010, Jung et al. 2011), or as recently shown, 
to a quantum valley Hall effect (Marino et al. 2015). Second, it does not consider the crystal 
symmetry of the lattice, which may give rise to crystalline topological insulators, and even 
more general behaviour (Slager et al. 2013). Last, but not least, it does not take into account 
out-of-equilibrium systems. 

The case of topological insulators (TI’s) under the influence of a time-periodic 
perturbation, the so-called Floquet TI’s (FTI’s) (Kitagawa et al. 2011, Inoue & Tanaka 2010, 
Gu et al. 2011, Lindner et al. 2011, Ezawa 2013, Fregoso et al. 2013, Rudner et al. 2013, Wang 
et al. 2013, Fregoso et al. 2014, Gomez-Leon et al. 2014, Perez-Piskunow et al. 2014, Quelle 
& Morais Smith 2014, Kundu et al. 2014, Usaj et al. 2014, Carpentier et al. 2015), has, 
until now, been considered in three unequal regimes. Firstly the so-called quasi-equilibrium 
regime, where J <C hu <C A; here J is the hopping parameter, which is roughly the 
bandwidth of the relevant set of bands, u> the driving frequency, and A the gap between 
the next set of bands. In this case, the system constituents cannot follow the perturbation, 
and the system remains at quasi-equilibrium with simply renormalised lattice parameters. It 
is the regime that has been most studied (Eckardt et al. 2005, Zenesini et al. 2009, Koghee 
et al. 2012, Struck et al. 2013, Goldman & Dalibard 2014). Secondly, the regime where J -C 
Hu ~ A. This regime is starting to attract interest in optical lattices (Parker et al. 2013, Zheng 
et al. 2014, Zheng & Zhai 2014), but it has been unexplored in the context of condensed matter. 
Thirdly, the regime where J ~ hu <C A, and where the equilibrium topological classification 
breaks down. It is this regime where most of the work on FTI’s in condensed matter has taken 
place, and these kind of systems have even been simulated in twisted photonic waveguides 
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(Rechtsman et al. 2013), where the third spatial dimension takes the role of time. There have 
been attempts to define Chern-type topological invariants valid for every frequency range 
(Lindner et al. 2011, Rudner et al. 2013, Carpentier et al. 2015), and it is known that these 
invariants reduce to the equilibrium ones in the first regime. The transition between the first 
and third regime has been investigated theoretically in Ref. (Kundu et al. 2014) for graphene 
irradiated by circularly polarised light. However, for the case of ultracold atoms, this regime 
has so far been overlooked. 

Here, we show that for circularly shaken honeycomb optical lattices the transition 
between the first and third regime can be understood in terms of band inversion. These 
band inversions occur because, with decreasing frequency, the different Floquet bands start 
to overlap. Because the circular shaking induces phonon resonances, avoided crossings that 
generally host topological edge states occur. The polarisation of the shaking breaks time- 
reversal symmetry, such that the resulting FTI is in the quantum anomalous Hall class. 
It has recently been shown that one-phonon resonances create an additional topological 
gap in the spectrum at non-zero energy, whereas two-phonon resonances destroy the 
topological nature of the zero-energy gap by creating counter-propagating edge states (Quelle 
& Morais Smith 2014). Here, we derive an effective continuum model from the out-of- 
equilibrium lattice model, and show that, in the vicinity of the band inversion occurring for the 
two-phonon resonances, it turns out to be described by the Bernevig-Hughes-Zhang (BHZ) 
model. Additionally, we demonstrate the validity of this model by comparing our results 
with a numerical solution of the full problem. In this way, our results describe the transition 
between the quasi-equilibrium and the resonant regimes of the shaken system, by modelling 
the appearance of phonon resonances in the latter regime in terms of band inversions. This 
provides a link between non-equilibrium and equilibrium topological states of matter. Finally, 
we discuss a possible experimental observation in shaken honeycomb optical lattices loaded 
with ultracold fermions, as well as in the recently realised honeycomb superlattices of CdSe 
nanocrystals (Kalesaki et al. 2014, Boneschanscher et al. 2014). 

2. The Hamiltonian in co-moving coordinates 

Consider an optical lattice that is shaken in time. The deviation of the lattice from its 
equilibrium position is denoted by r(t); by assumption r(t + T) = r(t ) for some period 
T. To find the Hamiltonian in co-moving coordinates, we consider the Poincare-Cartan form 

dS = p ■ dq — Hdt (1) 

along the trajectory of a system in phase space. We change coordinates to co-moving 
coordinates q = q + r(t), p = p, such that dp = dp and dq = dq + r(t)dt, where r(t) is the 
time derivative of r(t). The Poincare-Cartan form can thus be rewritten as 

dS = p ■ dq — Hdt = p ■ [dq — r(t)dt\ — Hdt (2) 

= p ■ dq — (H + p ■ r(t)) dt. (3) 

We immediately read off the Hamiltonian in co-moving coordinates: H = H + p ■ r(t). The 
extra term encodes the pseudoforces seen because the co-moving frame is not inertial. For a 
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shaken optical lattice, the Hamiltonian reads 



(4) 


where the potential determines the lattice, which we take to be honeycomb. 


In the co-moving frame q and p, Eq. (|4]) becomes 



(5) 


( 6 ) 


For circular shaking, r(t) = r 0 (cos(u;f),sin (cut)), which means ||r(f)|| 2 is constant in time 
and can be ignored by shifting the energy. The circular shaking thus induces a rotating vector 
potential eA, which has constant magnitude e/1 = mr 0 u. Compare this to the Hamiltonian 
for a system irradiated by circularly polarised light, where eA = eE/u for electric field E. 

3. Floquet theory 

The Hamiltonian in Eq. (|5]) is periodic in time, so according to Floquet theory (Sambe 1973, 
Hemmerich 2010), the time-dependent Schrodinger equation has quasi-periodic solutions 
0(f) = exp (—iet/h)(j)(t), where 0 is a periodic function in time and thus a solution of 
H F (j)(t ) = e0(f). Here, the Floquet Hamiltonian is defined as H F := H — ihd t . If H acts on 
the Hilbert space H, and Ht is the Hilbert space of T-periodic functions, then H F acts on 
H 8) Ht- The space H F is spanned by the functions | n) := exp(inut), and has inner product 
1 f dt/T. With respect to the states n), we can write H F as a block matrix with elements 



: H m _ n + rnhuoS 

m,n • 


(7) 


m—n 


Here, the H m _ n are the Fourier modes of the original Hamiltonian //. 

4. Model 

We apply the Floquet formalism to fermions in a circularly shaken honeycomb lattice, with 
shaking radius r 0 , and frequency oj. We work in the co-moving reference frame, where the 
Hamiltonian has the form in Eq. ([5]). To facilitate our analysis we use the tight-binding 
approximation, as done for graphene (Castro Neto et al. 2009). In second quantisation, 
the result is the Bloch Hamiltonian, except that one has to account for the vector potential 
according to the Peierls substitution, k k := k + eA/h : 



( 8 ) 


Here, J > 0 is the NN hopping amplitude, and our convention for the NN hopping vectors Si 
is <5 0 = a(0,1), and 6±i = —a (±a/ 3, l) /2, where a is the NN distance. Consequently, one 
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Hn=J 


(9) 


obtains, via Eq. (|7|), the matrices H n 

0 K -n ex p[* (k- Si + cq iH )] 

Kn exp[i (—k ■ Si + ai tn )\ 0 

In Eq. ([9]), K n := J n {aimypjf K ), where J n is the Bessel function of the n th kind, and 
ai !n := nArg[<^] + mr/2, where Arg gives the angle of a vector with the x-axis. Using 
Eq. (|9]), we obtain H F as an infinite block matrix 

( : : : : \ 

■ ■ ■ Hq hw 

H f = ••• H—i H 0 H x ■■■ . (10) 


H 1 

h 2 

Ho 

Hi 

H -1 

Ho — hw 


V ; 




If the spectrum around an energy nhw is desired, it can be computed by truncating the right- 
hand side of Eq. © around H 0 + nhw. By tuning the shaking frequency u, one may 
reach certain regimes where H F exhibits topological edge states. In the large-frequency 
limit, ftcc » J, the Floquet bands are well separated, but one finds a hierarchy of band 
crossings upon decrease of ui, when hw « ./. These band-crossings create additional, possibly 
topological gaps in the spectrum. In the following, we focus on the band crossing in the 
vicinity of k = 0 and e = 0, appearing in the interesting regime: hw « 2.9 J. In Fig. |Tj 
the spectrum of H F is plotted for hw = 3 J. The bottom of the valence band from H 0 + hw 
(on top) and the top of the conduction band from Ho — hw (below) are visible with a gap 
between them. As w is lowered, the valence band (on top) descends and the conduction 
band (on bottom) ascends; a band inversion takes place at hw ~ 2.9J, creating an avoided 
crossing [Fig. [2jA)]. At the avoided crossing, a single pair of edge states crosses the gap, as 
is highlighted in Fig. (2jb), where we show a zoom in on a narrow energy window. It should 
be noted that the edge states at k — 0, which correspond to the two-phonon resonance, are 
counter-propagating with respect to the zero-phonon resonance edge states at the Dirac points. 
(These also occur in the gap at s — 0, but have been made translucent to avoid confusion. 
They are depicted at e = ±hw, which are equivalent to e = 0 due to the periodicity of the 
Floquet spectrum; in these gaps, instead, the two-phonon resonance states have been made 
translucent.) The appearance of new edge states at k = 0 removes the topological protection 
of the edge states in the e = 0 gap (Quelle & Morais Smith 2014). 

This lack of topological protection can be better verified in a lattice with armchair 
termination, where the zero-phonon and two-phonon resonances both occur at k = 0. Indeed, 
they gap out because of hybridisation. Furthermore, if there is a domain wall in the system, 
where the orientation of the irradiation changes, the states localised at the domain wall also 
gap out due to hybridisation effects (Quelle et al. 2014). This shows that the appearance of the 
two-phonon resonance indeed makes the e = 0 gap trivial. The situation can be reversed by 
applying a staggered sublattice potential: this will destroy the zero-phonon resonance at the 
Dirac points, but leave the two-phonon resonance untouched. In this case, the appearance of 
the two-phonon resonance changes the gap at e = 0 from topologically trivial to non-trivial. 
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To write down an effective theory that allows us to characterise the band crossing, the gap and 
the dispersion of the edge states, we extract the relevant energy bands from H F . This is done 
by diagonalizing H 0 , and for simplicity we keep terms up to second order in k. We define the 
unitary transformation 


U ■= -y/2 ( 1 1 

2 ^ 1-1 

and consider the transformed Hamiltonian 

/ •• 


Hr = 


H 0 + hw H x H 2 
H-i H 0 //, 

H_ 2 H -1 H 0 - hw 


\ 


(ID 


\ 


( 12 ) 


/ 


where H n = UH n U by construction. From Eqs. (|9| and (11), one finds the identity 

3. 


Ho = Jkq ( 3^ -k ] o z + ... 


It follows that for hut ~ 3 J and km 0, the (2, 2) matrix element of the matrix H 0 + hu and 
the (1,1) element of H 0 — hu are much smaller than all the other energy scales in the problem. 
They are the zeroth-order energies of the two bands near the band crossing shown in Figs. |T] 
and[2| This leads us to define a second unitary transformation V that is characterised by the 



-n -n/2 0 n/2 it 

k 


Figure 1. (Colour online) The spectrum of the Floquet Hamiltonian Hp is shown. Plots were 
made for a ribbon geometry with zigzag edges, and k denotes the Bloch momentum along 
the length of the ribbon. Two periods of the spectrum of Hp are shown for huj = 3./ and 
mr 0 oj 2 a = J. The relevant feature is the impending gap closure at e = 0 and k = 0, when 
the Floquet bands n = 1 and n = — 1 overlap. To highlight this, we have made all bands, 
except for n = 1 and n = — 1, translucent. Red and blue represent different edges; the entire 
spectrum is spin degenerate. 
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(a) 

GJ 


gj/2 


-u/2 



Figure 2. (Colour online) (a) Same as in Fig. [T] but for hoj = 2.7 J. The relevant feature is 
the band inversion at e = 0 and k = 0, where the Floquet bands n = 1 and n = — 1 overlap, 
creating a gap with topologically protected edge states. To highlight this, we have made all 
bands, except for n = 1 and n = — 1, translucent, (b) A zoom in on the band crossing from 
(a) is provided to make the details of the gap visible. 


matrix elements (n\V\m) = a + 5 m>n - 1 + cr^d m . a+ \ and that permutes the basis vectors in such 
a way that the (1,1) elements of the 2x2 matrices along the diagonal of H F are interchanged 
with the (2,2) elements diagonally above. Here, a± = (a x ± ia y )/2. With respect to this 
basis, the Floquet Hamiltonian H F VII F V reads 


/ •• 


\ 


H f = 


H e ft + flu H \ H 2 

H -1 H eS //, 

H_ 2 H_ i H e g — Huj 


v : : : 


(13) 


where the Hi are more complicated matrices obtained by interchanging elements of the Hj, 
which do not need to be defined here. It follows that 


H cS = 


(H 0 + M 2 ’ 2 H 2 / \ 

Hbl (H 0 - hw) 1 ’ 1 J ’ 


(14) 


where the superscripts denote a specific element of the corresponding matrix. Using the 
method from Ref. (Goldman & Dalibard 2014), which works because H e s is smaller than huj, 
the corrections to this term can be calculated in terms of commutators of the Hi (for example, 
the first correction is [H _ x , //, \/uj). These terms are all of higher order in r 0 and 1/c o. For 
small u) and/or large shaking amplitudes r 0 , these higher-order terms become significant. As 
described in Ref. (Perez-Piskunow et al. 2015), which the describes the irradiated condensed 
matter analogue of this system, increasing r 0 can lead to phase transitions without additional 
band inversions, as the off-diagonal blocks become sizeable. For values of cu around the one 
at which the two-phonon resonance appears, however, this does not yet occur for the shaking 
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radii r 0 that we discuss. This can be seen from a comparison with numerical calculations, 
which shows that Eq. ( |~i~4] ) is sufficient to accurately describe both the gap size and the 
presence of the topological states; see Fig. |2jb) for example. 

Although our method can also be used to model the one-phonon resonance, the form 
of the effective Hamiltonian will be different, as can be seen from the different topological 
properties connected with this resonance. Higher-phonon resonances occur at low values of 
u and hence require the inclusion of higher-order terms in the effective Hamiltonian. Using 
the definition of H 2 , one obtains 


H 2 A = 77^2 (iky - k x ) + ..., 


and thus the BHZ Hamiltonian 

H cS = 

where 


M + Bk 2 a 2 A(k x — ik y )a 
A(k x + ik y )a —(M + Bk 2 a 2 ) 


M hu B 3 A 

~ = T~ 3k °’ 7 = 4 K °’ 7 = 


- Ko- 

2 


(15) 


(16) 


These expressions are correct up to order ( mr 0 ua/h ) 2 , and agree with the ones derived by 
Kundu et al. (Kundu et al. 2014) upon replacing eE/u by mr 0 u. The presence of the Bessel 
function of the second kind, J 2 , through n 2 in Eq. < [T6| ), shows that the opening of a gap at the 
band inversion is a second-order phonon process. It should be noted that for the NN-hopping 
J < 0, one finds that A in Eq. ( [T6] ) acquires an additional minus sign, but the spectrum 
remains unaffected. 


6. Edge states and gap size 


From the effective Hamiltonian in Eq. ( |T3| ), and following Ref. (Qi & Zhang 2010), one can 
derive an explicit solution for the edge state in the infinite half-plane. Using perturbation 
theory to linear order in k, the edge states then disperse as 


El = EAk = — Jka 

16 


mr 0 ua 

h 


+ 


i.e. the edge states have a velocity quadratic in the frequency u. 

From H c s, an expression for the gap size A can also be derived, and one obtains 


A 


3 hu nnr 0 ua\ 2 


+ ... 


(17) 


By substituting the parameter values hu = 2.7 J and mr 0 u 2 a = J into Eq. ( fl7| ) yields a gap 
size A = 0.033 J, which is in good agreement with the numerical results shown in Fig.|2jb). 


7. Conclusion 

In conclusion, we have investigated fermions in a circularly shaken honeycomb optical lattice 
in the intermediate regime, where the shaking frequency is on the order of the bandwidth 
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hoj « 3 J. In this particular regime, the system is characterised by a substantial overlap 
between the Floquet side bands, and a series of band inversions can be created that generally 
host topological edge states. We have concentrated on the crossing associated with two- 
phonon resonances, at hw « 2.9 J, and we have shown that the relevant effective continuum 
model is just the BHZ model for HgTe quantum wells. This allows for an understanding of 
the transition between the quasi-equilibrium regime and the resonant regime in terms of well- 
studied effective models, and especially in terms of band inversion, now between adjacent 
Floquet bands. 

Considering that the model Hamiltonian also describes a honeycomb lattice irradiated 
with circularly polarised light, the question remains whether the discussed effects can be 
observed in condensed matter. In this case, the phonon resonances become photon resonances, 
but the prior calculations remain valid, simply by replacing mr 0 uj by eE/co. A natural 
candidate would be graphene, but the relevant hopping parameter J = 2.8 eV and the 
NN bond length a = 1.4 A in graphene would require unphysically large frequencies 
beyond the THz regime, and a very high field strength of E ~ 5.3 • 10 10 V/m. A more 
promising candidate is a self-assembled honeycomb lattice of CdSe nanocrystals (Kalesaki 
et al. 2014, Boneschanscher et al. 2014), which hosts an s-band exhibiting a dispersion 
similar to that of graphene. The hopping parameter in these artificial structures depends on 
the diameter and the contact area of the nanocrystals. A hopping parameter J = 25 meV, 
that is roughly two orders of magnitude smaller than that in graphene, has been theoretically 
predicted for nanocrystals with a diameter of 3.4 nm (Kalesaki et al. 2014). By using light 
with E = 10' V/m and huj = 65 meV, a gap of 1.5 meV is obtained for these parameters, 
which is 6% of the hopping J. In Ref. (Wang et al. 2013), the Dirac states at the surface 
of a 3D topological insulator are irradiated by circularly polarised light, and the resulting 
photon resonance gaps are detected using ARPES. Although thermal excitation of the Dirac 
electrons is observed, it is possible to measure the Floquet spectrum before the states have 
been excited away from the bands. This, together with the predicted bandgap, implies that the 
two-photon resonance should be observable in the recently synthesised artificial superlattices 
of CdSe nanocrystals (Kalesaki et al. 2014, Boneschanscher et al. 2014), or in predicted 
similar structures (Beugeling et al. 2015). 

It would be much more natural to attempt to realise the Hamiltonian in Eq. ( [X5] ) through 
the use of optical lattices. Honeycomb lattices have been manufactured in the past (Soltan- 
Panahi et al. 2011), and they are more promising for two reasons. The first reason for this is the 
much larger lattice constant, compared to condensed matter systems: since the vector potential 
enters the Hamiltonian in the combination eAa , this allows for smaller vector potentials. The 
second reason is that the circular shaking described here creates a vector potential of the form 
eA = ?nr 0 u>, as opposed to E/u, so that increasing the frequency actually increases the 
vector potential, rather than supressing it. In graphene, for example, the required frequencies 
suppress the vector potential too strongly, resulting in the necessity of unphysically large 
electric fields. In contrast, taking a honeycomb optical lattice with NN hopping J and recoil 
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energy E r = h 2 k 2 /‘2rn, we can rewrite 

aeA ar 0 ujm 1 2 r 0 hu 

— = — = 2 ar °K k ~ = ' 2, ’7 l E r ’ 

To obtain the bandgap derived in the previous section would require shaking by a frequency 
hu = 2.7 J, at a radius r 0 = aE r /(140 J). For potassium atoms loaded in an optical lattice 
with wavelength k = 1064 nm, which corresponds to E r & 4410 Hz, the shaking would be 
at several kHz, with a radius of several tens of nm. This suggests that honeycomb optical 
lattices are a very promising candidate for realizing the topological states discussed here. The 
richness of shaking protocols, which are a hallmark of optical lattices, together with these 
encouraging results, promise that the up-and-coming field of atomtronics could be a prime 
candidate for the experimental investigation of Floquet topological insulators. 

The authors would like to thank Daniel Vanmaekelbergh and Michelle Burrello for use¬ 
ful discussions. The work by A.Q. and C.M.S. is part of the D-ITP consortium, a program 
of the Netherlands Organisation for Scientific Research (NWO) that is funded by the Dutch 
Ministry of Education, Culture and Science (OCW). 
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